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Abstract 



Two- and three-pion correlations are investigated in cases when disoriented 
chiral condensate (DCC) occurs. A chaoticity and weight factor are used as 
measures of two- and three-pion correlations, and the various models for DCC 
are investigated. Some models are found to yield the chaoticity and weight 
factor in a reasonable agreement with recent experimental data. 
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I. INTRODUCTION 



It is expected that chiral symmetry will be restored within the high temperature region 
generated by high-energy nuclear or proton collisions. The phase transition associated with 
the restoration may be of a nonequilibrium nature, caused by the rapid expansion and 
cooling of the high-temperature region. Through such a phase transition, the chiral order- 
parameter could have a direction in the isospin space, different from that in the true vacuum. 
This phenomenon is called disoriented chiral condensate (DCC) and has been intensely 
investigated in the last several years JIJ. The most prominent signature of DCC is expected 
to be that the fractional neutral-pion number obeys a specific probability distribution, a 
simple version of which is 

p(/)^, (i) 

derived by assuming that the direction of condensates is isotropic in the isospin space Jj],[l3| . 
Here, / is the fractional neutral-pion number, a ratio of the neutral-pion number and the 
total pion number, given by / = jj — , M n ° , M . The distribution predicts a large proba- 

iv-n-^ -\-l\ 7V _ -riV-TTQ 

bility that the number of the neutral pions is observed much less than that of the charged 
pions. Apparently such events, known as Centauro events @, are observed in cosmic-ray 
experiments. Accelerator experiments, such as the Minimax and WA98 experiments, have 
attempted to measure this distribution [0,0, but no clear signature of such a distribution 
has been observed. 

In relativistic heavy-ion collisions, two-pion correlations have been used to obtain infor- 
mation about the size and shape of the pion-emitting source, based on the Hanbury-Brown 
Twiss (HBT) effect. Three-pion correlations are expected to provide possibly new informa- 
tion about the source [||,[| . The first experimental data regarding the three-pion correlations 
has been recently reported J/J. In our previous paper ||, we investigated effects of a not 
completely chaotic source on the two- and three-pion correlations, without referring to any 
specific dynamics involved in the phase transition. If DCC should occur, a natural question 
would be how well effects of DCC could be observed in pion interferometry. In this paper, 
we address that question. 

It has been pointed out that the two-pion correlations yield anomalous values when 
the pions are emitted from a single coherent source exhibiting DCC P-jOU- Perhaps it is 
reasonable to assume the pion-emitting source associated with DCC to be coherent. As we 
will discuss in Sec. Ill, the pions emitted from such a source do not exhibit the HBT effect 
and do not serve as a means of identifying DCC. It would be naive, however, to assume 
that the pions would be emitted by a single coherent source exhibiting DCC and by nothing 
else. For example, there could be more than one coherent source associated with DCC, or 
coherent sources could appear with the background of a chaotic source. It is quite reasonable 
to expect that such a domain structure might occur at the appearance of DCC [11,11,1]. 



We investigate the effects on both the two- and three-pion correlations, of various possible 
characteristics of a source having domain structure. 

A note on terminology. As mentioned above, we treat a source or multiple sources of DCC 
as coherent. If a source is not associated with DCC, we will call it generic. A generic source 
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or generic sources can be either coherent or chaotic. A generic source that is combined with 
DCC sources will, however, be regarded as chaotic in this work. In our previous work ||, we 
investigated models in which coherent and chaotic generic sources are combined in various 
ways. The DCC models in Sec. IV have the same structure as the models in the previous 
work, except that a DCC source or DCC sources have been replaced by a coherent source 
or coherent sources, respectively. In figures, we will compare numerical results obtained by 
the use of the two groups of models. When we do so, we will introduce abbreviated names 
of the models for clear identification of curves in the figures. 

In Sec. |TJ, pion spectra and correlation functions are reviewed, and a chaoticity and 
weight factor are defined as measures of the strength of two- and three-pion correlations. 
We investigate correlations for a single DCC domain in Sec. |II, and for a more realistic 
pion source involving DCC in Sec. [TV[ In Sec.|V|, we compare our results with a recent 
experiment and present a summary of this work. 



II. PION SPECTRA AND CORRELATION FUNCTIONS 

In this section, for clarity, we define the pion spectra and correlation functions that will 
be used in the rest of this work. Our definitions are conventional. 
One-, two- and three-pion spectra are given by 

Wi(p) = (o^o*), (2a) 
^(PbP^Hi'i^O, (2b) 

wt(Pi,P2, P 3) = kMJ^MaA^ ( 2c ) 

where i, j, and k stand for +, — , or 0, specifying the pion charge such as 7r + , 7r_, or 7r , 
respectively. All momenta are on-shell, p° = Vp 2 + m 2 , above and hereafter. (...) is formally 
in terms of the quantum state, or Tr{p...} in terms of the density matrix, p. 
In this work we introduce various models in place of the formal definition. 
In terms of the spectra, we define correlation functions as 

Cn P uP2)= , (3a) 

C 3 J (pi,P2,P3) = —, j k ■ (3b) 
Wl(p 1 )W((p 2 )W^(p 3 ) 



As in [ 12 1 , in order to exclude the effect of multiplicity fluctuation, we use normalized 



correlation functions for identical pions: 



Cf(Pi,P2) = T rp^rCn PuP2 ), (4a) 
ct\puP2, P3 ) = ^^-2)) ^^^^' (4b) 



where 
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(hi) = J d 3 pWl(p), (5a) 
(hiihi - 1)> = J d 3 p 1 d 3 p 2 W 2 ll (p u p 2 ), (5b) 
(h t (h t - l)(hi - 2)> = J d 3 Pl d 3 p 2 d 3 p 3 Wt(PuP2,P3). (5c) 



Here, hi is the number operator, hi = J d 3 pa^a p . 

We define chaoticity, Xij(p), and weight factor, u>ijk(p), as 



\Ap) = {c 1 2^P2)-i) , (6) 

\ / Pl =p 2=p 



C l 3 h (PuP2,P3) ~ 1 - (C t 2(PuP2) + C J 2 k (p 2 ,p 3 ) + Cf(p 3 ,Pl) ~ 3) 

UijkVP) — — 



2VW(P1^) - l)(^f (P2,P3) " 1)(C"(P3,Pi) 



Pl=P2=P3=P 

(7) 



Ajj(p) and oJijkip) are measures of the strength of the two- and three-pion correlations. Note 
that Xij and Uijk are defined for zero relative momenta. Uijk has been extracted from the 
recent experiment for small relative momenta [0, but the difference is immaterial since it 
is a slowly varying function of the relative momenta []5]|| . When \j{p) is unity or zero, we 
call the pion-emitting source chaotic or coherent, respectively. For a chaotic source, u)ij k {p) 
is unity. 



III. A SINGLE COHERENT SOURCE ASSOCIATED WITH DCC 

In this section, we examine the case of pions that are emitted from a single coherent source 
associated with DCC. The correlation functions in this case can be derived in two different 
ways, by the isosinglet method and by the generating functional method. The isosinglet 
method has been previously used to obtain the anomalous values of the two-pion correlation 
functions ||13|| . The generating functional method Jl4| is a more systematic approach, and 



we will apply it to more complicated cases in Sec. IV. Though the two methods yield the 
same result in this simple case, we will show both for the sake of comparison. Both methods 
show that no HBT effect occurs in this case. 



A. The Isosinglet Method 

An isosinglet state vector of 2N pions is given by 
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m = , = J ' N \0), (8) 
<(2N+ 1)! 



where 

/ = 2a + a_ — <2q. (9) 
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I is an isosinglet operator |15|JT3| j9[1 . Since the creation and annihilation operators are of 



a single mode, the momentum variables are not explicitly shown here. The probability of 
finding 2n neutral pions is 



2n.2N 



(N\) 2 2 2N \2n)\ 
(2N + l)!(n!) 2 2 2 " 

1 (JV,n»l). (10) 



2^/N 



This large N and n limit is consistent with Eq. ([]]). The positive-pion spectra come out to 
be 

W+ = ~N, (11a) 
W+ + =^N(N-1), (lib) 
W 3 +++ = ^N(N-l)(N-2). (11c) 

In the large N limit, the two- and three-pion correlations for various combinations of pion 
charges defined in Eqs. (|3|) are then 

C++ = Cff= = | (12a) 

Cf = i (12b) 

C° 2 ° = ~, (12c) 

0^ = 0^ = 1, (12d) 

^3 ±±0 = ^3 ±T ° = ^ (12e) 
C 3 ±0 ° = |, (12f) 

C 3 00 = y, (12g) 

The anomalous values of the two-pion correlations, Eqs. ( |12a| )- (|12cj ), agree with those given 
previously 0. We see that the three-pion correlation functions also yield anomalous values. 

DCC gives correlation functions different from those of a coherent generic source. We 
note, however, that the normalized correlation functions, C's, are the same as those of a 
coherent generic source, all being unity. The reason for this is seen as follows: The isosinglet 
method is applied for a single mode (momentum), and yields the same W's for all modes. 
The normalization condition, Eq. (5), thus guarantees that all C's as defined in Eq. (4) 
are unity, the same value as those of the coherent generic source. The pions emitted from a 
single coherent DCC source do not then exhibit the HBT effect and do not serve as a means 
of identifying DCC. 
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B. The Generating Functional Method 



In order to achieve the isosinglet property of DCC, pion spectra should be averaged 
over all directions of isospin space |14j]. We evaluate correlation functions with c- number 
source-current formalism. 



In c-number source-current formalism, a generating functional method is useful [it]. We 



define the generating functional (slightly differently from [16]) as 

d 3 p 



G\zt(j,),zfy)\ 



(2tt) 3 • 2p° 



&*(P)J?(P) + *(P)4(P)) 



(13) 



where Ji(p) is the c-number source current for the pion charge, i = +, — , or 0. is 
the statistical average of the fluctuation of Ji{p)- Using this generating functional, we can 
obtain 



( a pi a P2 



U Pn U <?1 U <?2 



5 n+m G[z*{p), Zl {p)} 



<H(Pi) ■ ■ ■ ^z* in {p n )5z h {q x ) ■ ■ ■ 5z jm (q r , 



z;(p)=0 



Following |14]], we set the source current of each kind of pion, Ji(p), as 

MP) = J(p) n 3 

m =F m 2 



J±(p) = J(pV 



V2 



where n is a unit vector in isospin space. The ratio of neutral pions, /, is given by 



f = n 



.3- 



(14) 

(15) 
(16) 

(17) 



and \n±\ 2 = (1 — /)/2. rij is independent of momentum because we assume that all modes of 
momenta condense in the same direction in isospin space. Due to the isotropic distribution 
of the unit vector, the average that appears in Eq. ( |T3"D becomes 



(18) 



When we treat identical pions, the integrand of the right-hand side depends on only n| or 



1 77,-1- 1 . In such cases, we can replace the integral as 

/3, 



/^(W-l)- = iV(/>-> 



(19) 



where P(f) is defined in Eq. ([I]). 

The neutral-pion spectrum for DCC is obtained as 

W?(p)= f 1 P(f)dfn 2 3 W(p) 



\w{p\ 



(20) 



where 

w ® = vib 1 JW|2 - (21) 

For the charged pions, we obtain 

Wf(p) = W°(p). (22) 
Similarly, we can derive two-pion spectra, 

Wt ± ( Pl ,p 2 ) = W^( Pl , P2 ) = lw(p 1 )W(p 2 ), (23a) 

15 

W^( Pl , P2 ) = ±W( Pl )W(p 2 ), (23b) 
15 

W°°(p 1 ,p 2 ) = ±W(p 1 )W(p 2 ). (23c) 
Three-pion spectra are 

Wt ±± (Pi,P2,Ps) = W^( Pl ,p 2 ,p 3 ) = ^W( Pl )W(p 2 )W(p 3 ), (24a) 

W± ±0 ( Pl ,p 2 ,p 3 ) = W 3 ± ^°(p 1 ,p 2 ,p 3 ) = JLw( Pl )W{p 2 )W(p 3 ), (24b) 

W 3 ±00 ( Pl ,p 2l p 3 ) = ^W( Pl )W{p 2 )W(p 3 ), (24c) 

W^°( Pl ,p 2 ,p 3 ) = ^W(p x )W{p 2 )W(p 3 ). (24d) 

Equations (23) and (24) yield the same correlation functions as those obtained by the use 
of the isosinglet method, Eqs. ([i~2[ ). Furthermore, W(p) of Eq. (21) corresponds to the 
pion-number density for the mode (momentum) p, and W n 's are n-products of M^(p)'s. 
Accordingly, when the normalization condition of Eq. (5) is imposed, the Cs automatically 
become unity, which is the same value of C"s as in the coherent generic case. We thus see 
that though the correlation functions are momentum-dependent in this method while they 
are independent in the isosinglet method, both methods yield the same C"s and also C s = 1. 

When the pions are emitted from more than one coherent source with DCC, or they 
are also emitted from a chaotic generic source, the HBT effect appears and DCC could be 
identified. We will discuss these cases in the following section. 



IV. A SOURCE WITH DCC DOMAIN 

In relativistic heavy-ion collisions, it would not be realistic to consider all the pions 
being emitted from one large DCC source, as in the model in the previous section. In this 
section, we investigate three models that we expect to be more realistic. In the first model, 
the source consists of one DCC source and one chaotic generic source, while in the second 
model the source consists of multiple DCC domains. The third model is a combination 
of the previous two models, multiple DCC domains with one chaotic generic source. The 
HBT effect appears in these models and helps identifying the DCC signature in the pion 
interferometry. 
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A. A Partially Coherent DCC Source (PC-DCC) 

In this subsection, we introduce a partially coherent source with DCC [r3L|TT|. Source 
current, Ji(p), is separated into two parts as 

Mp) = Jd(p)^ + JI( P ), (25) 

where Jd{p) is the DCC source current and J l g {p) is the generic one. We assume that the 
DCC source is coherent and that the generic source is chaotic. The average in Eq. fll3|) 
becomes 

(■■■)j= J ^(m- 1 ) / ■■•nww.wwww. ( 26 ) 



i=i 



where Vi[Jg*(p), J l g {p)] is a distribution functional of J l g {p) and assumed to have a Gaussian 
form, as in Ref. |T6|] , so that the higher-order moment of Jg(p) is represented by the second- 
order moment, for example, 



(jT(pi)4(pi))j(jT(p2)4(p2))j + (JT(Pi)J i g (P2))j(Ji*( P 2)J i g (pi))j. (27) 



At first, we concentrate on the case of positive pions. We introduce here several real 
functions, /y, <f)y, and These functions are defined as follows [Q, 

fij exp(i^) = 1 = JUPi) J D{Pi) ( 28 ) 

(2tt)3^2p? • 2p° 

F± eX p(i$+.) ee ^ ^ / o Q Q (Jr(Pi)$(Pi))j}> ( 29 ) 

where 

(■ ■ = / ■ • j+(p))vj+(p)vj?(p). (so) 

Hereafter, for simplicity, we denote i^j" as Fy. We obtain the positive-pion spectrum, 

W+{p l ) = ^dfP{f)Wf{ Pl ) 

= \fn + F 11 , (31) 

where 

W}(pi) = h| 2 /ii + *ii- (32) 
The two-pion spectrum is given by 
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W+( Pl , P2 )= f P(f)df{w?( Pl )W+(p 2 )+R-}(p 2 ,p 3 )} 



where 



^/n/22 + W+( Pl )W+(p 2 ) + R + ( Pl ,p 2 ), (33) 



4(pi,p 2 ) = F{ 2 + 2\ ni \ 2 f 12 Fl 2 cos ($ 12 - 12 ) (34) 
R l ( Pl ,P2) = F{ 2 2 + j^ 12 F* 2 cos ($ 12 - 12 ) . (35) 



The first term in the right-hand side of Eq. (|33| ) does not appear in the case of the partially 
coherent generic source ||, signifying the new aspect in the partially coherent DCC case. 
The two-pion correlation function becomes 

c 2 ++(p„ P2 ) = 1 + i ( i - - + w ,f ( ;^) <36) 

where e(p) is the ratio of the generic pion number and the total pion number: 

<Pi) = lf F l F (37) 



In this case, the normalization of Eq. ( pEaj ) is not unity in the limit of the infinite relative 
momenta. This is because the term |(1 — e(pi))(l — e(p 2 )) is momentum-dependent. In the 
following, we proceed by neglecting that the momentum-dependence of e(p) is independent 
of p, so that the correlation function is normalized as in Eq. 



77++/ \ -1 . 1 R + (Pl,P2) (tiQ , 

ta.^ = 1+ 1+ . (1 ^ )2 „, 1+(pi) „, 1+(p2) - < 38 > 
The chaoticity for normalized correlation function, Eq. (f|), is 

This chaoticity is different from that of a partially coherent generic source: 

Xgeneric = e(2 - e), (40) 

which is evaluated in ||. In Fig. [I], we show the chaoticity as a function of the ratio, 
e, for a partially coherent DCC source (PC-DCC). For comparison, we also show in the 
figure the chaoticity for a partially coherent generic source (PC-G) that has been previously 
investigated g,|JTJ. 

The three-pion spectrum is evaluated as 
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w 3 +++ ( Pl ,p 2 , P3 )= f dfp(f)\w;( Pl )w+(p 2 )w+( P3 )+ J2 w;(p a )R+(p blPc ) 



(a,b,c) 



+ 2 



Fl 2 F 23 F 31 C0S($12 + $23 + $31 ) 

1-/ 



E fabF bc F ca cos((p ab + $ 6c + $ m ) I 
(a,6,c) ^ ^ 

^l + (Pl)^l + (P2)^l + (P3) + E ^l + (Pa)i? + (p6,Pc 

(a,6,c) 



+2 
1 



F 12 F 23 F 31 C0S($12 + $23 + $31 ; 



+ 3 E fabF bc F ca cos((f) ab + $ 6c + $„ 
(a,6,c) 

19 1 2 

+ TTjv/l 1/22/33 + 77 X/ faafbbF cc + 2J faafbcFbc COs($f, c 

»40 40 , , . 40 , , . 



?6c 



)• (41) 



The weight factor for normalized three-pion correlation functions is 

"3e(2-e) + 2e 2 (3-2 e ; 



a; 4 



1 + §(l- e )3 + |(l-6)^ 



fe(l-e) 2 , 

^ 



_3 



(42) 



Figure illustrates the weight factors for a partially coherent DCC source and a partially 
coherent generic source as functions of e. The weight factor of a partially coherent DCC 
model becomes negative for small e. This is the aspect of the DCC case that characteristically 
differs from the generic case: 



2e 



U), 



generic 



(2-e) 3 / 2 ' 



(43) 



which is always positive. 

The chaoticity and weight factor for normalized neutral-pion correlations are 



A 



00 



^000 



e(2 



l + |(l-6) 



3e(2-e) + 2e 2 (3-2e) + fe(l-e) 
l + f(l-6) 3 + f(l-e) 2 e 



3A 



00 



(44) 



(45) 



These are shown as functions of e in Figs. [I] and |2|. The chaoticity for the normalized 
correlations of differently charged pions vanishes, because the generic part exhibits no HBT 
effect and the DCC part is taken to be coherent. 

It should be noted that the distribution function of the fractional neutral-pion number 
in the partially coherent DCC model is different from Eq. (JJ). The fractional neutron-pion 
number in the partially coherent DCC model is 



w7°(pi) 



W/(p 1 ) + Wf(p 1 ) + W^ Pl 
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f-fn + Fn 
fu + SF n 

= /(I - e) + \e, (46) 

where each kind of pion emitted from the generic source is assumed to have the same 
multiplicity, or = = F{[. e is that of Eq. (|37|) . The new probability distribution 
function for f, P'(f') is obtained by the use of P(f)df = P'(f')df as 

' 'h < f < i - h) 



P'(f') = 2^{f'-\e){l-e) 3 3 . (47) 

y (otherwise) 

This distribution function becomes Eq. (fl|) when e = 0, or one DCC source with no chaotic 
source. When e = 1, or no DCC source, P'(f') is meaningful only at /' = |. The charge 
fluctuation caused by DCC is reduced as the generic source increases. 



B. Multiple DCC Domains (M-DCC) 

In actual experiments, we expect that DCC would appear not simply in a possibly large, 
unified area but in several, separated domains. In order to describe such a domain structure, 
we use multiple coherent sources whose number obeys the Poisson distribution. 



In an analogy to the description of a generic source in ||12|| , we write the source current 
for N DCC domains as 



N 

Mp) = T,j(Py P ' Xn ~ t(>nn n,i, (48) 
n=l 

where n n ^ is a unit- vector in isospin space, describing the direction of condensate in the n-th 
domain. The position of the n-th domain is X n , distributed according to the function, p(X n ). 
The normalization of p(x) is / p{x)dx = 1. 9 n is a random number uniformly distributed 
between and 2tt, so as to satisfy the chaotic property among the different domains. iV 
obeys the Poisson distribution, but the zero-domain event, or iV = 0, should be excluded 
from it. Instead of the usual Poisson distribution, the appropriately renormalized Poisson 
distribution is then 

N i 

(renorm) " 1 r Ar i / ac\\ 

lN iVTI^ foriV=l~oo. (49) 

In this case, the average of the generating functional, Eq. (0), is 

oo . / N j3 \ f / N 

(• • -h = e i { N norm) / n - 1) / n dx NP (x 

N=l J \n=l ^ / J \n=l 




(50) 



The pion spectra are obtained as 
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w; 



W+( Pl ) 



a \j{Pi)\ 2 

3 l _ e -a ' 

li(pi)l 2 li(p 2 )| 2 r2 



1 



15 



-a 



1 2 , -L 2 1 |2 

-a H — a pi9 
9 9 IP 1 



(51a) 
(51b) 



|j(Pi)| 2 |j(P2)| 2 |j(P3)P 
1 - e" a 
> 2 1 , 



35 
+ 



-a 



15 



-a 



27 



-a 



4 2 1 3 

— a H a d 

45 27 



EIp^I 2 + 



(a,6) 



27 



-a 



2Re(pi 2 p 2 3P3i) 



where 



Pij = / rf 4 xp(x)e- i(Pl ^ ) :c . 



(51c) 



(52) 



Using the above equations, we derive the chaoticity and weight factor for normalized corre- 
lations as 



A j 



a 



a + 



6 • 



« 2 + !« + fi 



'a + 



54 



a 



(53) 
(54) 



If the domains are not DCC sources but coherent generic sources, the chaoticity and weight 
factor become [01 



A 



a 



generic 



U). 



generic 



a + V 

12a 2 + 2a + 3 la + 1 
2 a 2 + 3a + 1 V a 



(55) 
(56) 



respectively. In Figs. |3] and f|, we show the chaoticity and weight factor of the multiple 
DCC model (M-DCC) as functions of the mean number of the DCC domains, (N) = a/(l — 
exp(— a)). For comparison, we also show in the figure the chaoticity and weight factor of 
the multiple coherent generic model (M-G) that we examined in our previous work ||. 

In this model, the correlations of differently charged pions has no HBT effect. For 
example, the two-pion correlation function of a paired positive pion and negative pion is 



W 2 + (pi,p 2 ) 



\j(Pl)\ 2 \j(P2)\ 



2 r 



— a H — a z 
15 9 



(57) 



W<2~~(pi,p2) is independent of pi2- Note that n 2 ^ appears during the calculation, but it is 
reduced to no contribution after averaging over the azimuthal angle by taking account of 
the isotropic distribution of condensate: 



1 i-2-k 

P(f)df / 
o Jo 



2tt 



-71, 



-^P{f)df J exp(2z0)g 



0. 



(58) 
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where n\ = y/1 — f cos0, n 2 = a/1 — / s i n 0- The HBT effect thus vanishes for a paired 
positive pion and negative pion by this averaging. Similarly, the correlation of a neutral 
pion and a charged pion shows no HBT effect. 

The chaoticity and weight factor of neutral pions are 

ct 

Aoo = (59) 



"°°°=2 a* + fa + f W - (60) 



These are illustrated as functions of the mean number of domains in Figs. ^ and 

As in the case of the partially coherent DCC of the previous subsection, the fractional 
pion number is also modified in the case of multiple DCC domains. We can write the 
probability distribution of / for iV DCC domains as 



PnU)= d[f--^ / f n \- 7 =---—=df 1 df 2 ..-df N , (61) 

where /„ is the ratio of the neutral pion number and the total pion number in the n-th 
domain. The pion spectra emitted from each domain is assumed to be the same. The 
observed distribution is obtained by averaging the above probability distribution with the 
Poisson distribution, Eq. fl49|) , 

a N 1 
jv=i iV • e 1 

Figure f| illustrates P a {f) for the mean number of the domains, (N) = 1,2,4, and 20. 
(N) = 1 corresponds to Eq. (|l]). Figure |5| is similar to a figure in where no detail of 
the derivation is given. The distribution indeed becomes sharper as (N) increases [ffTlffl. 



C. Multiple DCC Domains with a Chaotic Generic Source (M-DCC/G) 

We now consider the case of multiple DCC domains with one chaotic source emitting 
generic pions. This model corresponds to the combination of the two models in the previous 
subsections. Note that the Poisson distribution used in the current model is the usual one, 
instead of Eq. ©: 

^usual) = 0_ e „ a (g3) 

The source current is 

N 

Hp) = E i(p)e*- x -- tf -n B , i + 4(p), (64) 

n=l 

where the variables are defined in the same way as in the previous subsections. For conve- 
nience, we introduce S z ab and T a b. 
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rp J*(Pa)j(Pb) 
J-ab — — : l==Pab 



2(27r)Vpgpg 



(65) 
(66) 



where py is defined in Eq. ( |52|) and is equivalent to exp(i$^ fe ) in Eq. (p9|) . Hereafter, 
when we do not denote the superscript i explicitly, we imply that + is excluded. 
The one-pion spectrum and two- and three-pion correlations are 



Wt(pi) 

C+ + (p llP 2) 
Ct ++ {Pl,P2,P?>) 



a 



— Txi + Sii, 





3^12 + 5*12 


2 


h , 6(l-e)2" 
. ' 5 a 


W?{pi)W?{ V 2) 




(67a) 
(67b) 



1 i 54 (1-c) 3 , 18 (l-e) a 
~ r 35 a 2 5 a 

2 



f T ab + S ab + £a(l - e ) |TJ 2 + f(l - e) • 2Re(T a6 S ( 



6a J 



W?(p«)W?(p 6 ) 

(f T12 + 5*12) T23 + S23) ff T31 + S31 



Wf( P i)W+(p 2 )W+(p 2 ) 
where the ratio of the generic pion number and the total pion number, e, is 



(67c) 



S 
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-J- 11 



We assume that e is independent of p. The chaoticity and weight factor are obtained as 

a 



(68) 



A 



++ 



a + f(l-e) 2 



a 2 + |a(l - e) 2 + S(l - e) 3 (13 - 28e) , a + |(1 - e) 2 



a 2 + fa(l-e) 2 + |(l-^ 



51 , 



(69) 
(70) 



In Figs. |6| and |7|, the chaoticity and weight factor are shown as functions of the mean number 
of DCC domains, (N) = a, for various values of e. 

The chaoticity and weight factor of neutral pions are obtained as 



A 



a 



00 



^000 — 2" 



1 2a 2 + fa(l - e) 2 + - e) 3 (17 - 42e) a + |(1 - e) 3 



175 

a 2 + fa(l -e) 2 + f (1 - e) 3 



a 



(71) 
(72) 



where the spectra of neutral generic pions is the same as that of the positive ones, or 
Sii = S^. The correlation functions of nonidentical pions have no HBT effect, as in the two 
previous models. 
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If TV DCC domains are produced, the fractional neutral-pion number is actually 



1 NT U + S° n + S+ + S u 
_ En 3/ n (l-e) + ea 
3A(1 - e) + 3ea 

As a consequence, the probability distribution of / is 




(73) 



' e(/) = £^ ■/ V" 3A(1 - e) + 3ea: in^Tf ( 74 ) 



V. DISCUSSIONS AND SUMMARY 

As we have seen in Sec. Ill, the pions emitted from a single coherent source of DCC 
do not yield pion interferometry useful for identifying a signature of DCC. Also, when the 
pions are emitted from a source in which a single coherent DCC source and a chaotic source 
are combined, this situation remains practically the same. Figures [l] and |3] show that the 
partially coherent DCC source (PC-DCC) and the partially coherent generic source (PC-G) 
each yield a chaoticity and weight factor similar to that of the other. The difference between 
the two is numerically too close for interferometry to serve as a means identifying a DCC 
signature. If pion interferometry for the neutral pions should become feasible, combined 
data of charged and neutral-pion interferometry would serve as a means. Note that the pion 
interferometry of different charges exhibits no HBT effect and is not useful. 

Figure |8| illustrates the weight factor as a function of the chaoticity in the models exam- 
ined in this work, except for the model of multiple DCC domains with one generic source 
(M-DCC/G). For comparison, we also show in this figure the recent experimental data from 
the CERN NA44 Collaboration [0, the chaoticity and weight factor being 0.4 — 0.5 and 
0.20 ± 0.19, respectively. (We have combined the systematic and statistical uncertainties 
quadratically. ) The partially coherent DCC model (PC-DCC) is closest to the data point, 
while the partially coherent generic model (PC-G) is next. The ratio of the generic pion 
number and the total pion number, e, that yields the closest PC-DCC is found from Fig. 2 
to be about 0.30. That is, about 70% of pions are emitted from the DCC domains. This 
implies that we have a large charge fluctuation of the pions emitted. Such a large charge 
fluctuation does not seem to have been observed in the experiments H, however. Note that 
the essential aspect of the model that yields the results close to the data is the inclusion of 
a chaotic generic source. In Fig. |], we see the model of multiple DCC domains (M-DCC) 
with no chaotic generic source yielding the weight factor far away from the data point. 

Figure || shows that the multiple DCC domain model (M-DCC/G) successfully yields 
the chaoticity and weight factor in agreement with the data. The parameter values for the 
best fit are found to be a = (N) = 0.18 and e = 0.57 from Fig. 7. This implies that the 
mean number of the DCC domains is 0.18 and that the ratio of the generic pion number 
and the total pion number is 0.57. For these parameter values, the charge distribution, Eq. 
(0),is 
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P a , e {f) « 0.83P (/) + 0.15Pi(/). (75) 

Here, Po(f) is equal to 8(f — 1/3) in our simple model, while P%(f) is approximately the 
inverse square root of /, as in Eq. ([[]). Since P\(f) is suppressed by the factor of 0.15, 
Pa,e(f) of Eq- (|75p is a distribution dominated by a sharp peak with a slow-varying l/\/J- 
like background. In practice, the sharp peak can be replaced by a smoother function, such 
CIS db binomial distribution peaking at / = 1/3. 

We have thus shown that the NA44 data can be explained by the model of multiple DCC 
with a chaotic generic source. We have not demonstrated, however, that the NA44 data 
prove the appearance of DCC. Unfortunately, if we limit ourselves to pion interferometry, 
we require interferometry of the neutral pions in order to strengthen the case. For the above 
parameter values, the chaoticity and weight factor for the neutral pions are A o = 0.35 
and cuooo — —0.12, respectively. If interferometry for the neutral pions should be carried 
out, these values should signal the observation of DCC. Note that pion interferometry of 
differently charged pions does not serve for identifying DCC because the generic part exhibits 
no HBT effect and the multiple DCC part also has no HBT effect due to the azimuthal 
average in isospin space (as discussed in the case of the partially coherent DCC model in 
Subsec. |I V A| and in the case of multiple DCC domain model in Subsec. |I V B|) . 



In summary, we have investigated two- and three-pion correlations when DCC occurs. 
A chaoticity and weight factor are used as measures of two- and three-pion correlations and 
the various models for DCC have been investigated. The chaoticities and weight factors for 
DCC are different from the generic case. The existing experimental data are in agreement 
with the model of multiple DCC domains with one generic source. We suggest that the 
chaoticity and weight factor for neutral pions together with those for the charged pions 
should enable the identification occurrences of DCC. 
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FIGURES 




8 

FIG. 1. Chaoticity, A, as a function of the ratio of the generic pion number and the total pion 
number, e. The solid curve represents the partially coherent generic model (PC-G). The dashed 
and dotted curves represent the partially coherent DCC model (PC-DCC) for the positive pions 
and for the neutral pions, respectively. 
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FIG. 2. Weight factor, lo, as a function of the ratio of the generic pion number and the total 
pion number, e. The solid curve represents the partially coherent generic model (PC-G). The 
dashed and dotted curves represent the partially coherent DCC model (PC-DCC) for the positive 
pions and for the neutral pions, respectively. 




1 1.5 2 2.5 3 3.5 4 4.5 5 

<N> 

FIG. 3. Chaoticity, A, as a function of the mean number of domains, (N). The solid curve 
represents the multiple coherent generic sources (M-G). The dashed and dotted curves represent the 
multiple DCC domain model (M-DCC) for the positive pions and for the neutral pions, respectively. 
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FIG. 4. Weight factor, to, as a function of the mean number of domains, (N). The solid curve 
represents the multiple coherent generic sources (M-G). The dashed and dotted curves represent the 
multiple DCC domain model (M-DCC) for the positive pions and for the neutral pions, respectively. 
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FIG. 5. Probability distributions in the multiple DCC domain model, P a {f), as a function of 
the fractional neutral-pion number, /, the mean domain number, (N) = 1,2,4,20, (corresponding 
to a = 0, 1.59, 3.92, 20.0, respectively). 
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FIG. 6. Chaoticity, A, of the positive pions for the model of multiple DCC domains with one 
chaotic generic source, as a function of the mean number of DCC domains, (N). The curves are 
shown for the ratio of the generic pion number and the total pion number, e = 0.1 — 0.9, with step 
0.1 from down to up. 




<N> 

FIG. 7. Weight factor, u, of the positive pions for the model of multiple DCC domains with 
one chaotic generic source (M-DCC/G), as a function of the mean number of DCC domains, (N), 
varying the ratio of the generic pion number and the total pion number, e, from 0.1 to 0.9 with 
step 0.1, from down to up at (N) = 5. 
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FIG. 8. Weight factor as a function of chaoticity for the positive pions. The solid and dashed 
curves represent the partially coherent generic model (PC-G) and the partially coherent DCC 
model (PC-DCC). The dotted and dot-dashed curves represent the multiple coherent generic model 
(M-G) and the multiple DCC domain model (M-DCC), respectively. The data point is that from 
the CERN NA44 [7]. 
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FIG. 9. Weight factor, uj, as a function of chaoticity, A, for the positive pions, for the model of 
multiple DCC domains with one generic chaotic source (M-DCC/G), varying e from to 1. The 
lines from down to up correspond to the mean number of domains, (N) = 0.1, 0.2, 0.4, 0.6, 0.8, 
1.0, 1.5 and 2.0, respectively. The data is from the NA44 experiment [7]. 
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